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ABSTRACT 

We report the Fermi Large Area Telescope discovery of y-ray pulsations from the 22.7 ms pulsar A in the double 
pulsar system J0737— 3039A/B. This is the first mildly recycled millisecond pulsar (MSP) detected in the GeV 
domain. The 2.7 s companion object PSR J0737— 3039B is not detected in y rays. PSR J0737— 3039A is a faint 
y- ray emitter, so that its spectral properties are only weakly constrained; however, its measured efficiency is typical 
of other MSPs. The two peaks of the y- ray light curve are separated by roughly half a rotation and are well 
offset from the radio and X-ray emission, suggesting that the GeV radiation originates in a distinct part of the 
magnetosphere from the other types of emission. From the modeling of the radio and the y- ray emission profiles 
and the analysis of radio polarization data, we constrain the magnetic inclination a and the viewing angle f to be 
close to 90°, which is consistent with independent studies of the radio emission from PSR J0737— 3039A. A small 
misalignment angle between the pulsar’s spin axis and the system’s orbital axis is therefore favored, supporting 
the hypothesis that pulsar B was formed in a nearly symmetric supernova explosion as has been discussed in the 
literature already. 
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1. INTRODUCTION 

Pulsed GeV y- ray emission from more than 100 pulsars 9 10 
has been observed by the Large Area Telescope (LAT; Atwood 
et al. 2009) on board the Fermi satellite, launched in 2008 
June (see The Second Fermi Large Area Catalog of Gamma- 
Ray Pulsars; Fermi LAT Collaboration 2013, in preparation, 
hereafter 2PC). The current population of y- ray pulsars includes 
objects known from independent radio or X-ray observations 
and detected in y rays by folding the Fermi LAT data at the 
known rotational periods (e.g., Abdo et al. 2009a; Noutsos 
et al. 2011; Espinoza et al. 2013), and y-ray pulsars found 
through direct blind searches of the LAT data (Abdo et al. 
2009b; Pletsch et al. 2012) or radio searches for pulsars in 
unassociated y-ray sources (see Ray et al. 2012 and references 
therein). These pulsars are energetic (spin-down luminosities 
E = Atx 2 IP/P 3 > 10 33 erg s 1 , where P is the spin period, P 
is its first time derivative, and I denotes the moment of inertia, 
assumed to be 10 45 g cm 2 in this work) and are typically nearby. 
They therefore have large values of the “spin-down flux” E /d 2 
and of the heuristic detectability metric ~J~E /d 2 (Abdo et al. 
2010a). 

The double pulsar system PSR J0737— 3039A/B (Burgay 
et al. 2003; Lyne et al. 2004) consists of two radio-emitting 
neutron stars in a tight 2.4 hr orbit. Radio timing observations 
of the two pulsars provide high-precision tests of strong field 


9 Resident at Naval Research Laboratory, Washington, DC 20375, USA. 

10 A list of y-ray pulsar detections is available at 

https://confluence.slac. Stanford. edu/display/GLAMCOG/Public+List+of+ 
LAT-Detected+Gamma-Ray+Pulsars. 


gravity (Kramer et al. 2006). The 2.7 s pulsar J0737— 3039B has 
a low spin-down luminosity £ ~ 1.7x1 0 30 erg s 1 , three orders 
of magnitude smaller than that of the least energetic y-ray pulsar 
currently known, making this pulsar unlikely to be detectable 
by the LAT. On the other hand, the higher E of ~5.94 x 
10 33 erg s _1 — for which the Shklovskii effect (Shklovskii 1970) 
is negligible due to the system’s low transverse velocity — and 
the modest parallax distance of 1 150^ 22 6 o pc measured with very 
long baseline interferometry observations (Deller et al. 2009) 
make the 22.7 ms pulsar J0737— 3039A a credible candidate for 
a detection in y rays. 

Searching for high-energy pulsations from PSRs J0737— 
3039A and B allows exploration of regions of the P—P diagram 
that are devoid of y-ray pulsars (Abdo et al. 2010a): mildly 
recycled millisecond pulsars (MSPs) for the former (mildly 
recycled pulsars being pulsars that were only partially spun- 
up by accretion of matter from a companion star), and slowly 
rotating normal pulsars for the latter. Increasing the variety of 
y-ray pulsar types helps us understand the phenomenology 
of high-energy emission from these stars, and how emission 
properties evolve with, e.g., E or age (with the caveat that 
the evolution could be masked by effects arising from the 
geometry). Additionally, the modeling of pulse profiles at 
different wavelengths can yield constraints on the geometry of 
emission, in particular, the inclination angle of the magnetic axis 
(a) and the observer viewing angle (() with respect to the spin 
axis. Constraints on a and £ can also be determined by studying 
the radio polarization. Since the inclination of the binary system 
with respect to the line-of-sight is accurately known, 7 directly 
measures the misalignment between the pulsar’s spin axis and 
the orbital angular momentum, a quantity that gives important 
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clues about the formation of the binary system (for a discussion, 
see, e.g., Ferdman et al. 2008). 

We here report the discovery of y- ray pulsations from the 
22.7 ms pulsar J0737— 3039A, in 43 months of data recorded 
by the Fermi LAT. In Section 2 we describe the analysis of the 
LAT data and the resulting light curve and spectrum for PSR 
J0737— 3039A. In Sections 3.1 and 3.2, we model the offset 
radio and y-ray light curves, as well as the radio polarization, 
under various emission models, allowing us to constrain the 
pulsar’s geometrical characteristics. We conclude with a brief 
discussion on the implication of these results. 

2. LAT ANALYSIS 

We analyzed Fermi LAT events recorded between 2008 
August 4 and 2012 March 7 belonging to the “Source” class 
under the P7_V6 instrument response functions, and with zenith 
angles smaller than 100°. We rejected events recorded when 
the rocking angle of the telescope exceeded 52°, when the 
instrument was not operating in the science observations mode 
or when the data quality flag was not set as good. Analyses of 
the LAT data were carried out using the Fermi Science Tools 11 
(STs) v9r28p0. The data were phase-folded with the Fermi plug- 
in (Ray et al. 2011) distributed with theTEMPo2 pulsar timing 
package (Hobbs et al. 2006), and the ephemerides for PSRs 
J0737— 3039A and B published in Kramer et al. (2006). Since 
2008 March PSR J0737— 3039B has been invisible in the radio 
domain (Perera et al. 2010) because of the precession of its 
spin axis (Breton et al. 2008) causing its radio beam to miss 
the Earth, precluding timing measurements contemporaneous 
with the Fermi mission. Additionally, changes in the radio 
profile of PSR J0737— 3039B prevented Kramer et al. (2006) 
from constructing a coherent timing model for this pulsar. 
Our ability to accurately fold the Fermi LAT data using this 
ephemeris is therefore limited. Nevertheless, the very low spin- 
down luminosity of this pulsar makes it unlikely to be detected 
in y rays. 

From initial pulsation searches using standard data selection 
cuts we found no evidence of y-ray pulsations from PSR 
J0737— 3039B. However, selecting events found within 0°5 
from PSR J0737— 3039A and with energies above 0.2 GeV 
we obtained a value for the bin-independent FI- test parameter 
(de Jager & Biisching 2010) of 21.3, corresponding to a 
pulsation significance of ~3.7er. The y-ray pulse profile of 
PSR J0737— 3039A comprises two peaks, which is reminiscent 
of known y-ray pulsar light curves (e.g., Abdo et al. 2010a), 
with no evidence for emission in the [0; 0.2] U [0.5; 0.8] phase 
interval. 

As illustrated in Kerr (2011a) and Guillemot et al. (2012b), 
weighting the y-ray events by the probability that they originate 
from the putative pulsar increases the sensitivity to faint pul- 
sations. To measure these probabilities we performed a binned 
likelihood analysis of the spectra of the sources in the region of 
interest, using the PyLikelihood module of the STs. Our spec- 
tral model included the 50 2FGL catalog sources (Nolan et al. 
2012) within 20° of PSR J0737-3039A. The contribution from 
the pulsar, which is not associated with any 2FGL source, was 
modeled as an exponentially cutoff power law (ECPL) of the 
form dN/dE oc E r exp (— E / E r ), where T is the photon index 
and E c is the cutoff energy of the spectrum. The source model 
also included contributions from the extragalactic diffuse emis- 
sion and the residual instrumental background, jointly modeled 

11 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 


Table 1 

Main Properties and y-Ray Parameters of PSR J0737— 3039A 


Parameter 

Value 

Rotational period, P (ms) 

22.7 

Period derivative, P (10 18 ) 

1.76 

Spin-down luminosity, E (10 33 erg s -1 ) 

5.94 

Magnetic field at the light cylinder, #lc (10 3 G) 

4.97 

Distance, d (pc) 

1150+ 220 

1 1 JU _ 150 

First peak FWHM, FWHMi 

0.43 ± 0.01 

First peak FWHM, FWHMi 

0.03 ± 0.02 

Distance from the closest radio peak maximum, <5i 

0.20 ± 0.01 

Second peak position, O 2 

0.93 ± 0.01 

Second peak FWHM 2 

0 02 +0 03 
u ' u -0.02 

Distance from the closest radio peak maximum, <$2 

0.13 ±0.01 

y- ray peak separation, A = O 2 — O] 

0.49 ±0.01 

Photon index, T 

< 1.3 

Cutoff energy, E c (GeV) 

0.4 ± 0.4 (< 0.9) 

Photon flux, F (^0.1 GeV) (10 ~ 9 photons cm ~ 2 s _1 ) 

6 ± 3 (< 9) 

Energy flux, G (1^0.1 GeV) (10 -12 erg cm -2 s _1 ) 

4 ± 1 (< 5) 

Luminosity, L y /f, q = 47 T Gd 2 (10 32 erg s -1 ) 

67 3 2 (< 8 ) 

Efficiency, 17 //q = 4irGd 2 /E 

0.10!°om (< 0.13) 


Notes. See Section 2 for details on the measurement of these parameters. 
The quoted errors are statistical, while the numbers in parentheses represent 
canonical values, obtained by fixing the photon index of PSR J0737— 3039A at 
1.3 in the spectral analysis. 

using the iso _p7v6source template, and from the Galactic dif- 
fuse emission, modeled with the gal_2yearp7v6_v0 map cube. 12 
The parameters of PSR J0737— 3039A, of the seven sources 
within 10° of the pulsar, and the normalization factors of diffuse 
components were left free in the fit. A Test Statistic (TS; for a 
definition, see Nolan et al. 2012) value of 20.3 was found for 
PSR J0737— 3039A at this stage. 

To improve the quality of the spectral results, we analyzed 
the [0; 0.2] U [0.5; 0.8] phase interval (OFF pulse), removing 
the contribution from the pulsar. This allowed us to obtain an 
improved fit of the spectra of neighboring sources. Sources 
with TS values below 2 were excluded from the best-fit model 
obtained at this point, and the spectral parameters of sources 
beyond 3° from the double pulsar system were frozen at their 
best-fit values. We finally analyzed events in the complementary 
[0.2; 0.5] U [0.8; 1 ] phase interval (ON pulse) with the resulting 
source model, refitting the contributions of the pulsar and of the 
three sources within 3°, to determine the spectral parameters 
for PSR J0737— 3039A listed in Table 1. A TS value of 31.2 
is obtained for PSR J0737— 3039A. Refitting the data with a 
power law shape for the pulsar, we found that the ECPL model is 
preferred by the likelihood at the 4 o level. A cross-check of this 
analysis performed with the independent Pointlike analysis 
tool (Kerr 2011b) yielded results consistent with those listed 
in Table 1. Also listed in the table are the y- ray luminosity 
L y = 4jrfnGd 2 where G is the phase-averaged energy flux 
above 0.1 GeV and assuming a beaming correction factor 
fn= 1 , and the derived efficiency // = L y /E. The photon 
flux F and energy flux G of PSR J0737— 3039A are among the 
lowest values of any y-ray pulsars detected to date, probably 
a consequence of its low E value. We also note that PSR 
J0737— 3039A has the lowest value of the magnetic field at 
the light cylinder of any recycled pulsars detected in y rays, 
Blc ~ 4.97 x 10 3 G (see Equation (3.22) of Lorimer & Kramer 
2005). 

12 See http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 
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Figure 1. Radio and y - ray light curves for PSR J0737— 3039A. Bottom 
panel: 1.4 GHz radio profile based on observations made with the Green 
Bank Telescope. Upper panels: Fermi LAT profiles in different energy bands, 
obtained by weighting each event by its probability to have been emitted by PSR 
J0737— 3039A. Photons with probabilities smaller than 0.05 were excluded. The 
horizontal dashed lines indicate the estimated background levels (see Guillemot 
et al. 2012b for the determination of these lines). The y-ray pulse profiles have 
40 bins per rotation. Two rotations are shown for clarity. 

(A color version of this figure is available in the online journal.) 


None of the different attempts to model the spectrum of PSR 
J0737— 3039A presented above yielded satisfactory measure- 
ments of its photon index, 17 In all cases, a small value consistent 
with 0 was favored. This is likely a consequence of the cutoff 
energy E c being too low to reliably determine the photon index 
r. Limits on the other spectral parameters can thus be placed 
by repeating the analysis using a value for the photon index 
that is representative of the T values observed for other MSPs. 
Table 1 lists the results obtained by fixing T at 1.3, this value be- 
ing the average photon index measured for the strongest MSPs 
in 2PC. The spectral parameters quoted in parentheses represent 
canonical values, arising from the lack of constraints on T. 

The tool gtsrcprob was finally used to assign each event 
a probability that it originated from PSR J0737— 3039A based 
on the fluxes and spectra obtained from the likelihood analysis. 
The integrated pulse profile of PSR J0737— 3039A over 0. 1 GeV 
and for probabilities larger than 0.05 is shown in Figure 1. The 


zero of phase is defined by the maximum of the first Fourier 
harmonic of the signal transferred back to the time domain. We 
find a weighted H- test parameter of 41 (Kerr 201 la; de Jager & 
Busching 2010) corresponding to a significance of ~5.4a. 

As can be seen from Figure 1, the y-ray profile comprises 
two peaks, with indications for additional complexity. We fitted 
the profile above 0. 1 GeV with Lorentzian functions and found 
the peak positions and FWHMs for the first and second y- 
ray peaks as listed in Table 1. Defining <5j and 82 as the 
separation between the y-ray peaks and the maxima of the 
closest 1.4 GHz radio peaks, respectively at <t> r j ~ 0.24 and 
<D„2 ~ 0.79, we find <5j = 0.20±0.01 and <5 2 = 0.13±0.01. 
The uncertainty on <5i and 82 caused by the error on the dispersion 
delay is A(DM )/(kf 2 ), where A(DM) = 5 x 10~ 3 pc cm~ 3 is 
the uncertainty on the dispersion measure (DM) reported in 
Kramer et al. (2006), k is the dispersion constant (see Lorimer 
& Kramer 2005), and / = 1.4 GHz is found to be ~5 x 10 -4 in 
phase and can thus be neglected. The y- ray emission from PSR 
J0737— 3039A is thus offset from the radio emission, suggesting 
distinct origins in the magnetosphere of the pulsar. This also 
holds for the X-ray peaks detected by Chatterjee et al. (2007) in 
Chandra HRC data, which are likewise aligned with the radio 
peaks. 

Finally, we have searched the OFF-pulse fraction of the data 
for modulation at the orbital period, caused by the collision of the 
particle winds from the pulsars, and we have also searched the 
ON-pulse signal for attenuation of the emission from pulsar A 
around conjunction, caused by, e.g., photon-light-pseudoscalar- 
boson oscillation in the magnetosphere of B, as proposed by 
Dupays et al. (2005). In both cases, we observed only steady 
emission as a function of orbital phase. 

3. CONSTRAINING THE VIEWING 

GEOMETRY OF PSR J0737-3039A 

We can place constraints on the magnetic inclination a 
and viewing angle 7 by modeling the radio and y- ray light 
curves and the radio polarization. For all models we assume 
the vacuum, retarded-dipole magnetic field geometry (Deutsch 
1955). The emission is assumed to originate in the open zone, 
determined by field lines that do not close within the light 
cylinder (at Aj c = c/Q. = cP/ 2 tt), traced to foot points on the 
star surface at Rns (radius of the neutron star), which define the 
polar cap. The computation and fitting procedure are described 
in Watters et al. (2009), Venter et al. (2009), Romani & Watters 
(2010), and references therein. 

In two of the prevalent models of y- ray emission, uniformly 
emissive zones of the magnetosphere stand in for more physical 
models. The first is the outer gap picture (OG; Cheng et al. 1986; 
Yadigaroglu & Romani 1995) in which radiation is emitted 
between the null charge surface where St ■ B = 0 and R\ c. The 
second is the two-pole caustic model (TPC; Dyks & Rudak 
2003), which we take to be a geometric realization of the 
slot gap model (Muslimov & Harding 2004). In the original 
TPC model emission extends from the surface to 0.75 Aj,c- 
In both geometries the y- ray emission is confined toward the 
edge of the open zone coinciding with the last open field lines, 
in a “gap” idealized as a region of width w em , interior to an 
accelerating layer of width u; acc . In the OG model the y- ray 
emission originates from a thin layer on the inner edge of the 
gap, closest to the magnetic pole, whereas in the TPC model 
emission is produced throughout the gap. The radio emission can 
occupy a large fraction of the open field lines, albeit generally 
at lower altitude. 
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Figure 2. Top: modeled light curves and y - ray data for PSR J0737— 3039A. Bottom: 1.4 GHz radio profile and best-fitting profiles. TPC light curves are shown as 
pink lines and OG light curves are shown as green lines. The vertical dashed pink line and the dash-dotted green line indicate the closest approach to the magnetic axis 
under the best-fit TPC and OG models, respectively. 

(A color version of this figure is available in the online journal.) 


3.1. Joint Radio and y-Ray Light Curve Fits 

One approach is to fit the y -ray and radio light curves directly, 
following the geometrical light curve modeling of Venter et al. 
(2009), with the fitting procedure described in more detail in, 
e.g., Guillemot et al. (2012b). For these computations we assume 
P — 25 ms and P = lxl0~ 18 ss _1 , and treat the radio pulse 
as a classical single-altitude, hollow-cone beam centered on 
the magnetic axis as described by Story et al. (2007). For this 
analysis we start from a grid of models with 1° steps in both 
a and £, and fit the emission zone accelerating and emitting 
gap widths, w acc and w em , with steps of 2.5% of the polar cap 
opening angle (0pc ~ (G/^ns/c) 1 ' 2 ). For this computation we 
extend the emission to 0.95 Rye- 

The models are fit to an unweighted-counts y-ray light curve 
for PSR J0737— 3039A, constructed by selecting events from 
the data described in Section 2 found within 2° of the pulsar 
with energies above 0. 1 GeV and with probabilities of being 
associated with the pulsar larger than 0.05. In the case of low- 
statistics y-ray light curves with sharp peaks such as those 
observed for PSR J0737— 3039A, using ay 2 statistic with 
binned light curves may not be optimal as this statistic is often 
insensitive to the peaks. We have thus found that it is preferable 
to use the unweighted counts, which allow the use of Poisson 
likelihood, as opposed to the weighted counts which require 
using ay 2 statistic for binned light curves. The radio profile is 
fit with ay 2 statistic, assuming the same relative uncertainty for 
each phase bin, which is combined with the likelihood from the 
y-ray fit to determine the best-fit model parameters from a scan 
over the model phase space. We used 30 bins for both the y-ray 
and radio light curves in these fits. 

One difficulty in joint fitting of the radio and y-ray profiles 
is that the low statistical uncertainty of the radio light curve 


compared to the uncertainty of the y-ray bins causes the fit to 
be dominated by the radio data and to ignore the constraints 
provided by the y rays. Because our radio model is simplistic, 
this can lead to unrealistic solutions. We have therefore investi- 
gated different prescriptions for the “uncertainty” on the radio 
profile bins in order to balance the radio and y-ray contribu- 
tions to the joint likelihood. The following prescription gives 
the most satisfactory results. First, we select an on-peak interval 
for the y-ray light curve, which should dominate the likelihood 
and define o- gav e as the average, fractional uncertainty of the 
y-ray bins in this interval. We then define a radio uncertainty 
CT r = r max x cr gave , where r max denotes the maximum radio light 
curve bin. Although this prescription gives satisfactory results, 
it is arbitrary. We tested its sensitivity by varying er r by a factor 
of two (see Guillemot et al. 2012a); the resulting best-fit geome- 
tries changed by a < 4° (12°) and £ < 4° (9°) for the TPC (OG) 
model. The latter values can be considered systematic errors on 
our fit results. 

The best-fitting model light curves are shown in Figure 2, and 
the parameters are given in Table 2. The 1<7 uncertainties are 
estimated from two-dimensional likelihood profiles for a and 
C and one-dimensional profiles for the other parameters. The 
best-fit gap widths for both models are size 0, indicating that 
they are unresolved by our simulations. The simple geometric 
models we use neglect important physics, e.g., the effects of 
magnetospheric charges and currents on the gaps, and cannot 
perfectly fit the data. We use this discrepancy to set the scale of 
our parameter constraints by rescaling the likelihood by its best 
fit value, reduced by one-half the degrees of freedom (dof). In the 
Gaussian approximation, this is equivalent to setting the best- 
fit reduced y 2 = 1. With this normalization, the y 2 differences 
between the TPC and OG models are not significant. Also shown 
in Figure 2 are the closest approaches to the magnetic axis under 
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Table 2 

Radio and y-Ray Light Curve Fit Parameters for PSR J0737— 3039A 


Model 

a 


U>acc a 

^em 

— ln(likelihood) 

dof 

fa 


n 

(°) 

(% ©pc) 

(% ©pc) 




TPC 

?7 

O 

oo 

86 + _ 2 14 


0.0 ± 2.5 

121.0 

54 

0 - 893 & 

OG 

88-n 

742 ? 

0.0 ± 10.0 

0.0 ± 2.5 

123.1 

53 

0 - 953 $ 


Notes. Using these geometric models, the results are invariant under either a reflection across a = 90° or f = 90° . 
a In the TPC model the accelerating and emitting gaps are the same. 


the best-fit TPC and OG light curves. We note that our modeling 
finds different locations for the magnetic axis under the two 
geometries due to the ambiguity in defining the first and second 
y-ray peaks when the phase separation is close to 0.5. 

Using these geometric models, the light curves obtained when 
reflecting across a = 90° and £ = 90° are the same, but shifted 
by 0.5 in phase. We report confidence regions for a < 90° and 
£ < 90°, but note that identical, reflected regions exist in the 
other three quadrants. 

From the simulations we can estimate the beaming correction 
factor fa (Watters et al. 2009) which relates emission along a 
given line-of-sight to the total flux. The values for the models 
are given in Table 2 with estimated uncertainties based on the 
confidence regions in the a-t; plane. Both models suggest that 
this factor should be ~1, implying corrected y-ray efficiencies 
close to 10%. 

In addition to standard TPC and OG models, we have ex- 
plored two alternative geometries. In the first approach, a very 
low-altitude radio cone (low-altitude slot gap radio geometry) 
was invoked in conjunction with the usual TPC model for the 
y rays. This was done in the context of a radio cone producing 
peaks leading the caustic y-ray peaks. The second alternative, 
motivated by the idea that the radio may indeed have a dominat- 
ing leading peak with the radio profile lagging the y-ray light 
curve, assumed conal radio and y-ray geometries, with the y-ray 
cone being lower than that of the radio. The fits from these alter- 
native approaches were not satisfactory as they were unable to 
reproduce the observed radio and y-ray peak shapes and separa- 
tions, leading us to abandon these scenarios. However, they point 
to interesting avenues of model refinement, e.g., investigation 
of cones with non-uniform emissivities (such as patchy or one- 
sided radio cones), or non-aligned y-ray and radio cones. Such 
refinements, which are expected to lead to improved light curve 
fits, are beyond the scope of the current paper and will be devel- 
oped in future work (A. S. Seyffert et al. 2013, in preparation). 

We remark that the fitting results presented in this section 
were obtained by assuming that the observed radio emission 
from PSR J0737— 3039A originates from both magnetic poles. It 
has been proposed that the non- 180° separation of the two radio 
peaks and the overall symmetry of the profile imply that the 
radio emission originates from magnetic field lines associated 
with a single pole (see, e.g., Manchester et al. 2005). However, 
Ferdman et al. (2013) showed that this scenario would imply 
a high altitude origin for the radio emission in the pulsar’s 
magnetosphere. If both the radio and the y rays originate from 
the outer magnetosphere, we would expect the pulses to be 
aligned in phase (Venter et al. 2012). The two-pole origin for 
the radio emission therefore seems to be a more natural solution. 

3.2. Polarization Fitting 

When attempting to understand pulsar viewing geome- 
try, radio polarization can provide very powerful constraints. 


Although recycled pulsars are notoriously difficult to model, 
PSR J0737— 3039A shows substantial linear polarization struc- 
ture, and thus offers good prospects of constraining the orien- 
tation. In our polarization modeling we follow the conventional 
assumption that the electric vector position angle (P.A.) follows 
the projection onto the plane of the sky of the magnetic field 
line at the emission point. The P.A. sweep as one moves past the 
radio pole(s) thus probes the viewing geometry with particularly 
strong constraints on f — t, — a . Our modeling extends beyond 
the point dipole rotating vector model (RVM; Radhakrishnan & 
Cooke 1969) to follow the distortion of the magnetic field as the 
emission zone moves to a non-negligible fraction of R\ c (Craig 
& Romani 2012). In addition, we follow Karastergiou (2009) 
in treating “orthogonal mode jumps” for which the P.A. shifts 
by ±90°, and optionally account for the effects of interstellar 
scattering (negligible for PSR J0737— 3039A). For details, see 
H. A. Craig et al. (2013, in preparation). 

Figure 3 shows a polarimetric profile for PSR J0737— 3039A 
measured at 1.4 GHz with the Parkes radio telescope (R. N. 
Manchester 2012, private communication). The Stokes param- 
eters and Gaussian decomposition of the linear intensity com- 
ponents are displayed in the top panel, while the lower panel 
shows the P.A. data. In this section we refer to the brightest 
radio peak as P2. Both it and PI are multi-component. In po- 
larization fitting, these components are interpreted as different 
emission zones with different linear and circular polarization 
fractions and, possibly, different orthogonal mode states and 
emission altitudes. Progressively across the PI pulse we see a 
component with strong polarization at nearly constant P.A., a 
largely unpolarized peak, a component with rapid linear polar- 
ization sweep after the maximum and a separate, weakly linear 
polarized component on the flat pulse tail. For P2 the pattern ap- 
pears reversed, with a weak linear tail leading the pulse, a linear 
component with rapid P.A. sweep, an unpolarized peak and then 
a strongly linearly polarized component with little P.A. sweep. 
The flat P.A. components at the front of PI and the back of P2 
are not naturally produced by any component locked to a local 
magnetic field. We suspect that these components are controlled 
by magnetospheric plasma, either in fixing the P.A. directly or 
in strongly distorting the field at very high altitude. Here we fit 
the rapid P.A. sweep in the central linearly polarized component 
of each pulse. The parameters are a, f, fo (the phase of the 
total intensity peak with respect to the surface dipole magnetic 
axis), ^ 0 : the P.A. offset and the two emission heights r t and 
;'2 in units of If c- The fit values, the y 2 , and the full projected 
(multi-parameter) uncertainty ranges on the fit parameters are 
shown in Table 3. Figure 4 shows a projection of the y 2 surface 
of the fit onto the a — £ plane. 

The best-fit P.A. models for the two poles are displayed in 
the lower panel of Figure 3. Note that the weak trailing linear 
component in PI is separated from the central component by a 
zero in the linear intensity L, and so is very likely an orthogonal 
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Figure 3. Polarimetric profile for PSR J0737— 3039A as observed at 1.4 GHz with the Parkes radio telescope (R. N. Manchester, private communication). Top: Stokes 
parameter curves (black: I, red: L, blue: V) and Gaussian decomposition of the linear intensity components (dotted curve is all linear, the thick green curve gives the 
rapid sweep central component fit here). Bottom: position angle data (blue: central components, yellow: other PA. values, green: PI tail with an orthogonal mode 
jump). The smooth curves give the best fit model for the two poles while the red circles denote the boundaries of the open zone at the emission altitude. The arrows 
denote the phase of the closest approach of the magnetic axes to the Earth line-of-sight. The smallest impact parameter (5 is for PI at phase (J)b = 0 (phase of the 
magnetic dipole). 

(A color version of this figure is available in the online journal.) 


Table 3 

Polarization Fit Parameters for PSR J0737— 3039 A 


Otf 


h/Rlc 

ri/R-LC 

x 2 

dof 

A <p 

(°) 

(°) 






98 - 8 -1.5 

95 . 8 + 1 4 3 3 2 

0 01 +0 - 22 
u - Ui -o.oi 

0 11+0-49 
u - 11 -0.05 

48 

35 

0 443+0-008 
u.^H-j_oo55 


Notes. Errors are the extrema of the la contours in the full multidimensional 
parameter space. r\ is the emission altitude of the central component of PI and 
r 2 the altitude of the central component of P2. A</> marks the offset of the closest 
approach of the surface magnetic axis (in PI) from the total intensity peak in 
P2; this is phase 0.243 in Figure 1. 


mode jump. The jumped P.A. points cluster around the model 
prediction. We do not attempt to model the total intensity light 
curve or the P.A. behavior of the flat P.A. components leading PI 
and trailing P2 since these do not follow simple magneto spheric 
models. However, we can make some comments about the origin 
of this emission. First, the circles mark the phase ranges of the 
open zone at the altitude fit for the two central components. 
These components occupy the second half of the open zone 
at this altitude, terminating at the closed zone boundary, an 
excellent consistency check. However, the wings of the pulses, 
especially at large (p, lie beyond the circles and hence must either 
come from higher altitude or arise from the classical closed 
zone. Our picture of each pole is thus a cut through a hollow 
cone with the central component arising from low altitude and 
higher altitude components from the pulse wings. Since PI has 


48 54 60 66 


106 ° 

104 ° 
a 102° 

100 ° 

98 ° 

96 ° 

Figure 4. Projection of the y 2 surface of the RVM fit to the radio polarization 
data (central components) for PSR J0737— 3039A onto the a-f plane. The best 
fit orientation angles (o', f) are indicated by an arrow, while the contours are 
shown at the la, 2 a and 3er levels. 

(A color version of this figure is available in the online journal.) 



the smallest \fl\ — 3°, it is viewed at lower altitude than P2 with 
P = 14?6. 

These fits also determine the phase of closest approach of 
the surface dipole axis. This can be compared with predictions 
of the outer magnetosphere y-ray models at the fit a and £. 
We follow Romani & Watters (2010) in examining the a and 
C for the best fits to the y-ray light curve and then checking 
against polarization constraints. For the OG model the best 
light curve matches are found very close to the angles fit to the 
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70 


— 1 1 1 1 I 1 1 1 1 T 

J0737-3039A Light Curves 

at pol’n fit a = 99°, f=96° 


— OG 

- - TPC 

(both <50 = 0,12) 



0.5 1 1.5 

Magnetic Axis Phase 0 B 


Figure 5. Light curves for PSR J0737— 3039A at the a and f angles determined from the radio polarization study for the OG model (solid blue line) and the TPC 
model (dashed green line). The observed radio and y-ray profiles are shown as a solid red line and as a shaded gray histogram, respectively. See Figure 3 for the 
definition of the magnetic axis phase <f>g. The blue arrow denotes the location of the magnetic axis under the OG and TPC geometries. 

(A color version of this figure is available in the online journal.) 


polarization data. For the TPC model (here limited to 0.75 Rye), 
the best fits also occur near (90°, 90°). These are somewhat 
worse than the best OG values since they include too much 
unpulsed emission at all phases, but they are quite acceptable. 
In this regard the polarization fitting agrees well with the fits 
of Section 3.1, when transformed to the a > 90° and £ > 90° 
region. However, the very accurate phase measurements from 
the P.A. sweep highlight a tension between the two bands: the 
radio-determined pole is always ~0.1 in phase before the best- 
fit axis for the y pulse. This tension is also evident in the 
fits to the cruder pulse profiles of Section 3.1. However, non- 
vacuum effects provide a plausible explanation for this offset. 
In Kalapotharakos et al. (2012), it was noted that the primary 
effect on the high altitude pulse of increasing the magnetosphere 
conductivity was a growing lag in phase with respect to the 
surface dipole. Lags of S</> ~ 0.1 were seen for high conductivity 
magnetospheres. If we disconnect the y-ray and radio zero- 
phase definitions in the joint light curve fits, letting S be larger 
by as much as 0. 1 in phase, the best-fit a and f values changed by 
^4°, generally making the agreement between the two methods 
better. 

In Figure 5 we show the predicted OG and TPC y-ray light 
curves at the pulsar geometry determined from the polarization 
study, superimposed on the radio and the y- ray data, and using 
the phase of closest approach to the surface dipole axis as 
defined by the polarization fits. Both modeled light curves are 
acceptable, with the OG one being somewhat better (reduced 
y 2 of 1.45 compared to 1.56 for the TPC model). For these a 
and f angles we find beaming correction factors f(> of 0.92 for 
the OG model and 1.15 for the TPC model, similar to the ones 


given in Table 2 and therefore also implying y- ray efficiencies 
of -10%. 

4. DISCUSSION 

4.1. The Low Spectral Cutoff Energy and Magnetic Field at the 
Light Cylinder of PSR J0737—3039A 

In all outer-magnetospheric models, we expect the accelerat- 
ing electric field E\ \ to scale as (e.g., Abdo et al. 2010b): 

E\\ oc C(r)B LC wl cc , (1) 

with C(r ) some function of radius, and w acc the accelerating 
gap width. On the other hand, if leptons are accelerated in the 
radiation-reaction-limited regime, i.e., when the gain in energy 
due to acceleration by the electric field is balanced by curvature 
radiation losses, we expect that the spectral cutoff energy should 
scale as: 

E c oc E\(% X J 2 ex Cir^w^B^py 2 , (2) 

with p c the radius of curvature (at the position where the photons 
of maximal energy are produced). If we take C(r), p c and w acc 
constant for all MSPs (e.g., independent of P and P), this leads 
us to expect: 

E c oc Blc- ( 3 ) 

One does see a trend when plotting published values of E c 
versus B L q for all known gamma-ray MSPs, but it is much 
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weaker than expected, and the index is much smaller than 0.75, 
implying that our assumption of constant CAr), p c and w a cc 
across the MSP population is too simplistic. PSR J0737— 3039A 
is a mildly recycled MSP, having a relatively long P and hence 
the lowest Z?lc a / j3/2 / l/2 of all currently known gamma- 
ray MSPs. Assuming that the above empirical trend may be 
extrapolated down to low B\ q values, one would expect a 
relatively low cutoff energy. This is indeed seen, although PSR 
J0737— 3039A’s cutoff energy falls significantly below even the 
extrapolated trend. 

Apart from invoking the weak empirical trend discussed 
above, one may more generally observe that the relatively 
long period of PSR J0737— 3039A implies a lower acceleration 
potential and hence a lower spectral cutoff than for typical MSPs, 
as is observed from the spectral analysis. 

4.2. Constraints on the System’s Formation History 

From the modeling of the radio and y-ray light curves of 
PSR J0737— 3039A under the OG and the TPC geometries, 
and the modeling of radio polarization data (see Section 3), we 
constrained the magnetic inclination angle a and the viewing 
angle £ to be close to 90°. The best-fit a and £ angles obtained 
from the polarization study (see Table 3) agree well with 
those obtained from the light curve modeling (Table 2). The 
conclusion that a and f are close to 90° is in line with the 
results of Ferdman et al. (2013), based on the analysis of six 
years of radio observations of PSR J0737— 3039A, revealing 
no significant variations of its radio profile with time. Ferdman 
et al. (2013) find a and £ values of 90i2+ 1 1 6 6 3 2 and 9(T8+° 0 2 4 7 6 , 
respectively. These values were calculated assuming that radio 
emission is seen from both magnetic poles, that the orbital 
inclination is 88°7, and averaging the results from all radio 
pulse heights used (30% to 50%) and all observation epochs. 

By virtue of the relationship between £ and the misalignment 
angle between the spin axis and the orbital angular momentum, 
Sso (cf. Equation (3.36a) of Damour & Taylor 1992), our values 
for C close to 90° imply small Sso values unless the pulsar is at 
a special precession phase. Nevertheless, the fact that Ferdman 
et al. (2013) find a 95% confidence upper limit on Sso of ~3?2 
from their observations taken between 2005 June and 2011 June 
indicates that the pulsar was very unlikely to be at a special 
precession phase during the Fermi observation. Our constraints 
on £ therefore provide an independent confirmation that Sso 
is close to 0°. Since the spin axis of pulsar A and the angular 
momentum of the binary are likely to have become aligned 
through the accretion of matter from pulsar B’s progenitor star 
by pulsar A, a value close to 0° for the misalignment angle Sso 
supports the scenario under which a small kick was imparted to 
the system by the supernova explosion of pulsar B’s progenitor 
star. An electron-capture supernova, resulting from a collapsed 
O-Ne-Mg core as proposed by Podsiadlowski et al. (2005) could 
provide such a formation scenario for pulsar B (see Ferdman 
et al. 2013 and references therein for a detailed discussion). 

Because PSR J0737— 3039A is a faint y-ray emitter, we were 
not able to investigate any evolution of the y-ray profile as 
a function of time caused by a variation of the line-of-sight 
angle, £, in the four years of LAT data considered here. Since 
the misalignment angle Sso h as been shown to be close to 0°, 
changes in £ are likely small but are still expected. An exciting 
prospect, when significantly more Fermi LAT data are available, 
will be to monitor the evolution of the y-ray profile of PSR 
J0737— 3039A. To date, no y-ray pulsar light curves have even 
been observed to vary with time. 
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